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1.  INTRODUCTION 


The  current  technique,  recognized  worldwide  as  the  best 
available,  for  manually  estimating  tropical  cyclone  intensity 
based  solely  upon  satellite  imagery  was  developed  by  Vernon 
Dvorak  (1975).  Dvorak  recognized  the  relationship  of  the 
cyclone's  banding  features  (BF),  and  the  areal  extent  and 
configuration  of  the  central  dense  overcast  (CDO)  with  the 
intensity  of  the  storm.  The  technique  gives  a  subjective, 
although  systematic,  estimate  of  the  current  intensity  and  the 
expected  intensity  change  which  is  based  upon  a  "T-number."  The 
T-number  which  may  range  from  T1  to  T8,  is  an  estimate  of,  where 
the  cyclone  "fits,"  in  Dvorak's  model  of  the  chronological  life 
cycle  for  "average"  cyclones.  The  current  intensity  (Cl)  number 
is  the  same  as  the  T-number  when  the  cyclone  closely  follows  the 
modeled  amount  of  daily  cloud  feature  change.  Dvorak  also  shows 
a  relationship  between  the  Cl  number,  wind  speed  and  the  minimum 
sea  level  pressure. 

Upon  examination  of  many  hurricanes  near  the  southeastern 
coast  of  the  United  States  Senn  and  Hiser  (1957)  found  that  in 
general  the  spiral  rainbands  observed  on  radar  closely  resemble 
logarithmic  spirals.  They  developed  a  method  of  determining  the 
storm  center  using  a  logarithmic  spiral  overlay  fitted  to  the 
spiral  bands  of  the  radar  imagery.  The  accuracy  of  this  method 
is  comparable  to  that  of  conventional  synoptic  methods.  Hughes 
(1952)  combined  data  from  a  number  of  reconnaissance  flights  into 
large  Pacific  cyclones  to  obtain  a  generalized  pattern  of  winds 
at  low-levels  (—1000  ft).  He  showed  that  the  trajectories  of  the 
low-level  inflow  were  logarithmic  spirals,  and  computed  averaged 
profiles  of  wind  speed  (both  radial  and  tangential),  relative 
vorticity  and  divergence.  Lahiri  (1981),  under  the  assumption 
that  the  cloud  spirals  are  coincident  with  the  streamlines  of  the 
wind  field,  shows  that  the  spiral  constant,  a,  in  the  logarithmic 
spiral  equation,  r  =  EXP[a9],  is  equal  to  the  negative  of  the 
ratio  of  the  radial  and  tangential  velocities.  He  also  shows 
that  the  ratio  of  divergence  and  relative  vorticity  are  also 
equal  to  the  negative  of  the  spiral  constant. 


Infrared  satellite  image  counts,  from  which  equivalent 
blackbody  temperatures  of  cloud  tops  are  calculated,  have  been 
shown  to  contain  meaningful  information  concerning  Intensity  and 
intensity  changes  associated  with  tropical  cyclones  (Gentry  et 
al.,  1978).  Latent  heat  released  when  air  with  high  equivalent 
potential  temperature  ascends  into  the  deep  convective  cells  of 
tropical  cyclones  provides  the  primary  energy  source  for  the 
storms.  The  amount  of  latent  heat  released  is  related  to  the 
depth  of  the  convective  cells  and  their  associated  colder  cloud 
tops,  which  ultimately  is  converted  to  kinetic  energy  by  means  of 
a  complex  nonlinear  scale  Interaction  conversion  process.  The 
time  Involved  for  this  conversion  from  latent  heat  energy  to 
available  potential  energy  and  then  to  kinetic  energy  is  used  to 
explain  the  observed  time  lag  between  the  maximum  convective 
activity  and  the  maximum  storm  Intensity. 

The  results  of  the.  previous  studies,  itemized  above,  suggest 
a  physical  relationship  between  satellite-observed  cloud  patterns 
associated  with  tropical  cyclones  and  dynamic  parameters  known  to 
be  related  to  tropical  cyclone  intensity  characteristics.  As 
noted  by  tahiri  no  direct  method  exists  for  objectively  estima¬ 
ting  the  spiral  cloud  band  parameters  from  satellite  images. 

The  Spiral  Linearization  Technique  (SLT)  was  devised  to 
objectively  analyze  digital  satellite  image  data  so  that  cloud 
bands,  which  typically  approximate  spiral  shapes,  could  be  more 
directly  related  to  tropical  cyclone  characteristics,  such  as 
intensity.  Estimates  of  tropical  cyclone  intensity,  intensity 
changes,  movement,  and  position  are  currently  determined 
subjectively  from  satellite  images  and  other  forecast  aids.  The 
accuracy  of  these  estimates  is  highly  subjective  and  dependent  on 
the  analyst's  experience.  The  development  of  objective  techni¬ 
ques  to  analyze  and  extract  information  from  digital  satellite 
data  of  tropical  cyclones  is  essential. 


2.  THE  CURRENT  SIT 


The  Spiral  Linearization  Technique  software  was  designed  to 
utilize  the  Satellite-Data  Processing  and  Display  System's 
(SPADS)  unique  capabilities  of  allowing  interactive  digital  image 
processing  (Schramm,  Zeleny,  Nagle,  and  Weinstein,  1982).  The 
SLT  consists  of  the  following  three  main  steps.  The  flow  diagram 
for  the  SLT  is  shown  in  Figure  1. 

2.1  IMAGE  PREPARATION  ( SATMG) 

Upon  selecting  and  displaying  the  desired  tropical  cyclone 
image,  the  user  must  enhance  the  image  so  that  the  dominant 
spiral  feature  is  isolated.  The  center  of  the  spiral  feature  is 
then  defined  by  positioning  a  cursor  on  the  spiral  center  by 
means  of  the  graph  pen/pad  or  console  entry.  Then  subsetting  by 
means  of  an  inner  and  outer  radii  is  performed  next.  The  inner 
radius  is  chosen  so  that  the  CDO  is  encompassed  to  avoid  unneces¬ 
sary  pixel  transformations  in  step  2.  The  outer  radius  is  set  to 
include  the  entire  spiral  band.  This  is  shown  in  Figure  2. 

2.2  TRANSFORMATION  TO  SPIRAL  COORDINATES  ( SLTMG) 

The  preprocessed  image  can  be  simultaneously  transformed  to 
any  or  all  of  the  four  spiral  types  which  include  archimedian, 
hyperbolic,  logarithmic  and  parabolic  spirals.  The  transformed 
image(s)  is  displayed  and  examined  to  see  if  the  proper  zero 
degree  reference  angle  is  acceptable.  The  zero  degree  reference 
angle  is  the  starting  theta  ( e)  for  the  transformation  from  x,y 
space  to  r,e  space.  The  above  scheme  is  illustrated  in  Figure  3. 
The  transformation  algorithm  has  been  modified  to  allow  for  a 
variable  offset  radius,  R0,  defined  by  the  user,  which  allows  for 
a  more  realistic  spiral  pattern  configuration  due  to  the  fact 
that  spiral  cloud  bands  appear  to  be  spiraling  into  the  eyewall 
(of  radius  Ro)  rather  than  into  an  absolute  center  point.  The 
effect  of  an  offset  radius  in  the  spiral  configuration  is  shown 
in  Figure  4. 
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Figure  2.  Subsetting  by  inner  and  outer  radii 
before  transformation  to  spiral  coordinates. 
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2.3  PROCESSING  THE  TRANSFORMED  IMAGE  ( X4MG) 


The  transformed  image  must  be  processed  further  prior  to 
entering  the  statistics  routine.  All  pixels  of  the  Image 
contained  within  the  Inner  and  outer  radii  have  been  transformed, 
leaving  cloud  features  not  related  to  the  spiral  band  remaining 
In  the  Image.  The  Image  must  be  further  subsetted  so  that  only 
the  pixels  constituting  the  spiral  band  are  retained.  Once  this 
is  accomplished  the  statistical  evaluation  of  the  spiral  band  can 
begin.  The  “goodness  of  fit"  of  the  observed  spiral  band  to  the 
given  spiral  type  Is  measured  by  the  linearity  of  the  transformed 
spiral  band.  The  statistics  routine  computes  a  linear  least 
squares  line,  correlation  coefficient,  sum  of  the  errors,  root 
mean  square  error,  mean  error  as  well  as  the  coefficient  of 
variation.  The  SLT  slope  of  the  linear  least  squares  line  Is 
then  correlated  with  the  known  spiral  constants  of  several 
perfect  logarithmic  spiral  images  generated  with  the  SLTVY 
program,  and  is  found  to  obey  a  power  law  formula.  The  equation 
is  of  the  form,  Y  *  AXB,  where  Y  is  the  spiral  constant,  X  is 
the  SLT  least  squares  slope,  and  A  and  B  are  constants  with 
values  of  4.5838  and  -0.99214,  respectively.  This  relationship 
is  illustrated  in  Figure  5.  The  logarithmic  spiral  constant  can 
be  shown  to  be  related  to  the  inflow  angle,  a,  with  the 
relationship. 


TAN  a  =  f  To  =  a 


dr 

de 


(1) 


The  relationship  between  the  inflow  angle  and  the  SLT  least 
squares  slope  is  shown  in  Figure  6.  Logarithmic  spirals  with 
offset  radii  obey  a  slightly  modified  inflow  angle  relationship, 
i . e. , 


TAN  a 


1  dr  _  V 

•  r  d0  '  a'1  '  r'* 


(2) 


From  this  equation  It  is  seen  that  for  radii  near  the  offset 
radius,  the  inflow  angle  approaches  zero,  while  for  large  radii, 
r  >>  R0,  the  Inflow  angle  approaches  that  of  a  spiral  with  offset 
radius  equal  to  zero 
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The  calculation  of  a  one-dimensional  Fast-Four ler-Transf orm 
(FFT)  of  the  average  theta  value  at  each  radius  (in  pixel  units) 
along  the  least  squares  line  is  the  next  SLT  option.  A  graphical 
plot  of  the  power  spectra  values  (the  sum  of  the  squared  real  and 
Imaginary  FFT  coefficients)  vs.  wave  number  can  be  constructed  if 
desired.  Inferences  on  cloud  spacing  and  orientation  can  be  made 
by  examining  the  power  spectra  coefficients,  provided  the 
transformed  spiral  band  Is  free  from  geometric  distortion. 

The  final  processing  option  is  the  calculation  of  a  two- 
dimensional  FFT  of  the  entire  spiral  band.  The  SLT  has  the 
capability  of  displaying  the  magnitude  of  the  FFT  coefficients  in 
image  format.  To  compensate  for  the  rapid  decrease  of  the  Image 
spectra  with  increasing  frequencies,  the  displaying  of  the 
natural  log  of  one  plus  the  magnitude  of  the  FFT  coefficients, 

1  .e. , 

D(u,v)  *  loge  (1  +  | F( u  , v)  | ) ,  (3 ) 

greatly  facilitates  the  visual  interpretation  of  the  two- 
dimensional  FFT  display  (Gonzalez  and  Wintz,  1977)  as  illustrated 
in  Figure  7.  A  plot  of  the  radially  averaged  profile  of  the  two- 
dimensional  FFT  power  spectra  coefficients  was  also  devised.  The 
DC  component  is  used  as  the  center.  This  will  give  a  linear 
representation  of  the  behavior  of  the  mean  wave  numbers  in  the 
two-dimensional  FFT  result. 
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3.  DATA  AND  ANALYSIS 

The  dependent  data  set  consisted  of  two  mile  infrared  satel¬ 
lite  images  from  the  geosynchronous  GOES  WEST  satellite. 

Cyclonic  storm  intensity  estimates  and  positioning  were  given  by 
Gunther  (1980)  using  the  Dvorak  technique  of  satellite  cyclone 
analysis  (Dvorak,  1975).  This  data  set  was  supplemented  by  Fleet 
Numerical  Oceanography  Center's  archived  tropical  cyclone  data 
set  which  contain  maximum  wind  speed  at  six  hour  Intervals  for 
the  life  cycle  of  each  storm. 

Two  tropical  cyclone  cases  were  selected  from  the  limited  IR 
image  sets,  and  upon  the  availability  of  continuous  wind  and 
image  data.  The  two  cases  chosen  were  eastern  Pacific  hurricanes 
Dolores  and  Enrique.  Descriptive  life  histories  of  the  storms 
will  follow. 

3.1  HURRICANE  DOLORES,  17-23  JULY,  1979 

The  sixth  tropical  disturbance  of  the  1979  season  originated 
approximately  350  n  mi  south  of  the  Guatemalan  coast  and  began  a 
westward  progression  at  about  10  kt  on  14  July.  Intensifying 
cyclonic  circulation  by  0600  GMT  17  July  was  evident  and  the 
disturbance  was  upgraded  to  a  tropical  depression.  The  rate  of 
intensification  increased  and  by  1800  GMT  the  same  day  the 
depression  was  upgraded  to  tropical  storm  Dolores  near  11.2°N  and 
107. 0°W.  Sea  surface  temperatures  of  29°C  were  conducive  to 
further  rapid  development  as  the  storm  moved  around  the  southwest 
edge  of  an  upper  level  anticyclone  centered  over  Baja, 

California.  Wind  speeds  increased  to  65  kt  by  1800  GMT  18  July. 
The  initial  satellite  image  was  for  1345  GMT  18  July  and  was  the 
first  of  twenty-one  images  analyzed,  the  last  being  at  2145  GMT 
22  July.  Wind  speeds  continued  to  increase  and  reached  their 
maximum  intensity  of  105  kt  by  0000  GMT  on  21  July.  Dolores 
continued  to  move  northwest  and  began  to  weaken  over  cooler 
(25°C)  water  and  was  downgraded  to  a  tropical  storm  by  0000  GMT 
22  July  near  20.3°N,  121. 8°W.  By  0000  GMT  23  July  winds  had 
decreased  to  less  than  30  kt.  Figure  8  shows  Dolores  at  maximum 
intensity 
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8.  Hurricane  Dolores  near  maximum 
kt)  at  0215  GMT  21  July  1979. 
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3.2  HURRICANE  ENRIQUE,  17-24  AUGUST  19 79 


Nearly  a  month  elapsed  before  the  next  tropical  disturbance 
began  significant  development  near  11°N,  107°W.  Sea  surface 
temperatures  of  30°C  favored  continued  development  as  the 
disturbance  moved  westward.  At  1800  on  the  17th  the  disturbance 
was  upgraded  to  a  tropical  depression  and  six  hours  later  was 
upgraded  to  tropical  storm  Enrique  near  11.2°N,  114. 9°W.  Enrique 
reached  hurricane  strength  by  0000  GMT  19  August  near  12.8°N, 

119. 0°W,  but  then  weakened  as  the  storm  moved  over  27°C  water. 
Enrique  was  downgraded  to  a  tropical  storm  with  winds  of  55  kt  on 
the  21st.  By  1700  GMT  on  the  21st  Enrique  again  reached  hurri¬ 
cane  strength  and  began  a  rapid  intensification.  By  0000  GMT  on 
the  22nd  wind  speeds  reached  110  kt  and  reached  maximum  intensity 
of  125  kt  by  1200  GMT.  Satellite  images  beginning  at  1415  GMT  21 
August  and  ending  at  0415  GMT  23  August,  for  a  total  of  nineteen 
images,  were  analyzed  via  the  SLT.  Moving  northwestward,  Enrique 
began  to  move  over  cooler  water  and  was  downgraded  to  a  tropical 
storm  by  1200  GMT  23  August  near  20.4°N,  132°W.  Enrique's  final 
position  was  21.3°N,  135. 8°W  at  1800  GMT  on  the  24th  when  little 
more  than  a  low-level  vortex  remained.  Figure  9  shows  Hurricane 
Enrique  near  maximum  intensity. 

3.3  THE  CASE  STUDIES 

While  carrying  out  the  case  studies,  it  became  apparent  that 
the  logarithmic  spiral  type  best  conformed  to  the  spiral  configu¬ 
rations  evident  in  the  satellite  data,  and,  in  order  to  save 
computation  time,  was  the  only  spiral  type  employed.  If  a 
combined  cylindrical  vortex  and  sink  are  considered,  then  the 
stream  function  is 

i|»=-me  +  <loger,  (4) 

where  m/tc  is  the  spiral  constant  for  the  logarithmic  spirals 
traced  out  by  constant  values  of  ip  (Mi  1  ne-Thompson ,  1960). 
However,  the  selection  of  the  logarithmic  spiral  type  pre-empted: 
(a)  applying  the  one-dimensional  FFT  to  the  subsetted  transformed 


Figure  9.  Hurricane  Enrique  near  maximum  intensity 
(125  kt)  at  1415  GMT  22  July  1979. 
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spiral  bands  and  (b)  deducing  the  spectrum  of  cloud  sizes  due 
to  the  distortion  of  the  transformed  Image  caused  by  the  nature 
of  the  transformation  algorithm. 

The  approach  during  the  case  studies  was  to  process  each  of 
the  Images  with  the  SLT  In  an  Identical  manner.  While  some 
subjectivity  was  entered  into  the  analysis  by  the  selection  of 
enhancement  limits,  band  to  evaluate,  etc.,  the  storms  were 
treated  Identically  after  the  preprocessing  of  the  image  was 
accomplished.  To  keep  minor  statistical  variations  to  a  minimum, 
an  arbitrary  inner  radius  of  30  pixels  was  chosen  based  upon 
preliminary  experiments. 

The  SLT  parameters  which  were  recorded  for  each  image  are: 
maximum  IR  count  in  storm,  maximum  IR  count  In  the  spiral  band, 
minimum  IR  count  in  band,  average  IR  count  in  band,  and  the 
spiral  inflow  angle.  The  minimum  IR  count  in  the  band  was  the 
most  subjectively  defined  parameter  as  it  was  chosen  by  incre¬ 
menting  the  lower  enhancement  limit  until  the  spiral  band  became 
"prominent."  The  average  IR  count  in  band  was  defined  as  the 
mean  value  of  the  maximum  and  minimum  values  rather  than  an 
actual  mean  of  the  pixel  counts  comprising  the  band.  The  spiral 
Inflow  angle  recorded  was  that  for  large  radii  (R>>  R0). 

The  maximum  wind  speed  data  was  plotted  vs.  time  and  the 
wind  speed  values  at  the  times  of  the  satellite  images  were 
interpolated  from  these  plots.  A  three  point  running  mean  was 
applied  to  the  SLT  parameters  to  smooth  some  of  the  variance  out 
of  the  data.  SLT  data  for  both  Enrique  and  Dolores  vs.  maximum 
wind  speed  is  shown  in  Table  1.  The  rate  of  intensity  change  was 
calculated  using  a  forward  difference  scheme  with  the  user 
supplying  the  Initial  maximum  wind  speed  and  subsequently  using 
the  above  maximum  wind  speeds.  The  time  increment  was  the  time 
between  satellite  images.  The  rate  of  intensity  change  vs.  the 
SLT  data  for  Dolores  and  Enrique  are  shown  in  Table  2. 


Table  1.  SLT-derlved  parameters  for  Dolores  and  Enrique  vs. 
maximum  Mind  speed. 

DOLORES 


IRMAX 

IBMAX 

IBAV 

IBMN 

a 

ALPHA 

VMAX 

228.0 

228.0 

218.0 

208.0 

.575 

36.3 

44.000 

230.0 

222.0 

210.0 

197.3 

.641 

32.6 

48.000 

230.7 

218.3 

204.5 

190.0 

.610 

31.4 

66.000 

230.3 

213.3 

198.2 

182.3 

.639 

32.5 

66.000 

226.0 

215.3 

201.8 

188.3 

.638 

32.5 

66.000 

225.3 

214.3 

203.0 

192.3 

.598 

30.7 

68.000 

223.7 

213.0 

202.7 

192.3 

.558 

29.0 

70.000 

224.7 

212.0 

201.8 

191.7 

.581 

30.0 

80.000 

222.7 

203.7 

191.3 

179.0 

.517 

27.0 

94.000 

223.7 

205.3 

192.0 

178.7 

.588 

29.8 

98.000 

223.7 

204.0 

186.3 

168.7 

.505 

26.1 

100.000 

223.7 

210.3 

190.8 

174.7 

.647 

32.3 

102.000 

223.0 

206.7 

187.5 

169.7 

.501 

26.0 

104.000 

221.7 

204.3 

185.8 

167.7 

.606 

30.6 

98.000 

220.0 

194.0 

176.0 

155.0 

.509 

26.6 

96.000 

218.0 

196.0 

176.7 

156.0 

.675 

33.6 

92.000 

218.0 

200.0 

181.0 

161.7 

.583 

29.7 

90.000 

217.0 

208.0 

191.7 

175.0 

.643 

32.2 

88.000 

216.7 

200.7 

189.8 

172.7 

.569 

29.5 

84.000 

215.3 

205.7 

166.7 

167.3 

.552 

28.6 

76.000 

216.0 

206.5 

184.0 

161.0 

.504 

26.5 

70.000 

ENRIQUE 


223.0 

191.0 

181.0 

171.0 

.678 

31.8 

56.000 

222.7 

189.7 

180.5 

171.3 

.641 

32.4 

58.000 

222.7 

194.3 

186.2 

178.0 

.688 

34.6 

60.000 

223.0 

202.7 

193.8 

185.0 

.778 

38.2 

62.000 

222.7 

206.8 

199.0 

192.0 

.700 

35.1 

66.000 

222.3 

207.7 

199.5 

191.3 

.665 

33.3 

78.000 

222.0 

206.0 

197.2 

188.3 

.700 

34.5 

83.000 

221.3 

207.7 

192.8 

178.0 

.683 

33.8 

96.000 

221.3 

204.7 

184.8 

165.0 

.579 

29. r 

105.000 

219.7 

203.3 

183.3 

163.3 

.409 

22.3 

112.000 

218.7 

198.7 

180.5 

162.3 

.451 

24.0 

124.000 

216.7 

202.7 

185.5 

168.3 

.588 

30.0 

125.000 

216.3 

200.7 

183.2 

165.7 

.759 

36.9 

125.000 

215.7 

199.0 

180.5 

162.0 

.734 

35.9 

125.000 

215.0 

189.7 

172.7 

155.7 

.664 

33.1 

125.000 

213.3 

192.0 

171.3 

150.7 

.562 

29.3 

125.000 

210.3 

1914.0 

172.5 

151.0 

.589 

30.4 

124  000 

209.7 

200.7 

177.8 

155.0 

.699 

34.6 

1 1  7.  ’00 

209.0 

200.5 

178.3 

156.0 

.757 

36.8 

104.1  ;o 

All  SLT 

parameters 

are  3  pt. 

running 

means  VMAX 

(kt). 

ALPHA  (deg) 

Table 

2.  SLT - 

derived  parameters  for 

Dolores 

and  Enrique 

vs . 

rate 

of  intensity  change 

• 

DOLORES 

IRMAX 

IBMAX 

IBAV 

I BMN 

a 

ALPHA 

RIC 

228.0 

228.0 

218.0 

208.0 

.575 

36.3 

2.000 

230.0 

222.0 

210.0 

197.3 

.641 

32.6 

1.667 

230.7 

218.3 

204.5 

190.0 

.610 

31.4 

1.071 

230.3 

213.3 

198.2 

182.3 

.639 

32.5 

.000 

226.0 

215.3 

201.8 

188.3 

.638 

32.5 

.000 

225.3 

215.3 

203.8 

192.3 

.598 

30.7 

.556 

223.7 

213.0 

202.7 

192.3 

.558 

29.0 

.833 

224.7 

212.0 

201.8 

191.7 

.581 

30.0 

2.083 

222.7 

203.7 

191.3 

179.0 

.517 

27.0 

2.917 

223.7 

205.3 

192.0 

178.7 

.588 

29.8 

1.667 

223.7 

204.0 

186.3 

168.7 

.505 

26.1 

.417 

223.7 

210.3 

190.8 

174.7 

.647 

32.3 

.833 

223.0 

206.7 

187.5 

169.7 

.501 

26.0 

.417 

221.7 

204.3 

185.8 

167.7 

.606 

30.6 

-1.250 

220.0 

194.0 

176.0 

155.0 

.509 

26.6 

-.833 

218.0 

196.0 

176.7 

156.0 

.675 

33.6 

-1.667 

218.0 

200.0 

181.0 

161.7 

.583 

29.7 

-.833 

217.0 

208.0 

191.7 

175.0 

.643 

32.2 

-.417 

216.7 

206.7 

189.8 

172.7 

.569 

29.5 

-1.667 

215.3 

205.7 

186.7 

167.3 

.552 

28.6 

-2.222 

216.0 

206.5 

184.0 

161.0 

.504 

26.5 

-5.000 

ENRIQUE 

223.0 

191.0 

181.0 

171.0 

.670 

31.8 

.833 

222.7 

189.7 

180.5 

171.3 

.641 

32.4 

.833 

222.7 

194.3 

186.2 

178.0 

.688 

34.6 

.833 

223.0 

202.7 

193.8 

185.0 

.778 

38.2 

.833 

222.7 

206.0 

199.0 

192.0 

.700 

35.1 

1.667 

222.3 

207.7 

199.5 

191.3 

.665 

33.3 

5.000 

222.0 

206.0 

197.2 

188.3 

.700 

34.5 

2.083 

221.3 

207.7 

192.8 

178.0 

.683 

33.8 

5.417 

221.3 

204.7 

184.8 

165.0 

.579 

29.1 

3.750 

219.7 

203.3 

183.3 

163.3 

.409 

22.3 

2.917 

218.7 

198.7 

180.5 

162.3 

.451 

24.0 

2.500 

216.7 

202.7 

185.5 

168.3 

.588 

30.0 

.417 

216.3 

200.7 

183.2 

165.7 

.759 

36.9 

.000 

215.7 

199.0 

180.5 

162.0 

.734 

35.9 

.000 

215.0 

189.7 

172.7 

155.7 

.664 

33.1 

.000 

213.3 

192.0 

171.3 

150.7 

.562 

29.3 

.000 

210.3 

194.0 

172.5 

151.0 

.589 

30.4 

-.278 

209.7 

200.7 

177.8 

155.0 

.699 

34.6 

-2.917 

209.0 

200.5 

178.3 

156.0 

.757 

36.8 

-5.417 

3  pt.  running  means  RIC  (kt/Hr) 
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3.4  DEVELOPMENT  OF  PREDICTIVE  EQUATIONS 

The  equations  developed  rely  solely  upon  IR  satellite 
information  to  predict  maximum  wind  speeds  and  rate  of  intensity 
changes.  A  stepwise  multiple  regression  was  performed,  and 
depending  upon  the  significance  level  for  deleting  variables, 
retained  either  three  or  four  independent  variables.  The  three 
variables  retained  were:  maximum  IR-storm  (IRMAX),  maximum  IR- 
band  (IBMAX),  and  the  spiral  inflow  angle  (ALPHA).  The  fourth 
independent  variable,  when  used,  was  the  average  IR-band  (IBAV). 
Separate  regression  equations  were  developed  for  both  Dolores  and 
Enrique. 

The  equations  developed  to  predict  the  maximum  wind  speeds 
work  well  with  three  independent  variables.  The  addition  of  the 
fourth  predictor  only  accounts  for  an  additional  0%  of  the 
explained  variance,  in  the  case  of  Enrique,  while  for  Dolores 
there  is  virtually  no  change  in  the  percent  variance  explained. 
The  analysis  of  variance  tables  for  both  the  three  and  four 
independent  variable  cases  are  shown  in  Table  3  for  Dolores  and 
Table  4  for  Enrique.  The  regression  equation  for  Dolores  using 
three  independent  variables  is, 

Vmax  =  BO  +  B1*IRMAX  +  B2*IBMAX  +  B3*ALPHA  (5) 

where  BO  =  352.4,  B1  =  0.6148,  B2  =  1.869,  and  B3  =  -0.6001  are 
the  linear  regression  coefficients.  The  four  variable  equation 
is. 


Vmax  =  BO  +  B1*IRMAX  +  B2*IBMAX  +  B3*IBAV  +  B4*ALPHA  (6) 

where  BO  =  322  3,  B1  =  0.6816,  B2  =  -1.543,  B3  =  -0.2702  and 
B4  =  -0.6188.  Plots  of  the  predicted  vs.  observed  maximum  wind 
speeds  via  the  three  and  four  variable  regression  equations  for 
Dolores  are  shown  in  Figures  10  and  11. 
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Table  3.  Analysis  of  variance  table  for  Dolores  (maximum  wind  speed) 


THREE  INDEPENDENT  VARIABLES 

Source  of 
Variation 

Regression 

Residual 

Total 

**  GOODNESS  OF  FIT  ESTIMATES  ** 

Percentage  of  Variance  Explained  =  .67633E  02 
Standard  Deviation  of  Residuals  »  .10806E  02 
As  a  Percentage  of  Response  Mean  =  .13348E  02 


Degrees  of  Sums  of  Mean  P  (Exceeding 

Freedom  Squares  Squares  F-Value  F  under  HO) 

.30000E  01  .41479E  04  .13826E  04  .11841E  02  .19819E-03 

•17000E  02  .19850E  04  .11677E03 

. 20000E  02  . 61330E  04 


Regression  .40000E  01 
Residual  .16000E  02 
Total  . 20000E  02 


FOUR  INDEPENDENT  VARIABLES 
.41578E  05  .10394E  04 

.19752E  04  .12345E  03 

.61330E  04 


.84199E  02  .74351E-03 


**  GOODNESS  OF  FIT  ESTIMATES  ** 

Percentage  of  Variance  Explained  =  .67794E  02 
Standard  Deviation  of  Residuals  =  .11111E  02 
As  a  Percentage  of  Response  Mean  =  .13725E  02 


Table  4.  Analysis  of  variance  table  for  Enrique  (maximum  wind  speed) 


THREE  INDEPENDENT  VARIABLES 


Source  of 
Variation 

Degrees  of 
Freedom 

Sums  of 
Squares 

Mean 

Squares 

F-Value 

Regression 

. 30000E  01 

.94710E  04 

.31570E  04 

.12075E  02 

Residual 

.15000E  02 

•39217E  04 

•26145E  03 

Total 

.18000E  02 

.13393E  05 

**  GOODNESS  OF  FIT  ESTIMATES  ** 


Percentage  of  Variance  Explained 
Standard  Deviation  of  Residuals 
As  a  Percentage  of  Response  Mean 


.70718E  02 
.16169E  02 
.16429E  02 


Regression  .40000E  01 
Residual  .14000E  02 
Total  .18000E  02 


FOUR  INDEPENDENT  VARIABLES 
.1 0238E  05  .2 5534E  04  .11358E  02 

.31549E  04  .22535E  03 

.13393E  05 


**  GOODNESS  OF  FIT  ESTIMATES  ** 


Percentage  of  Variance  Explained  =  .76443E  02 
Standard  Deviation  of  Residuals  =  .15012E  02 
As  a  Percentage  of  Response  Mean  =  . 15253E  02 


P  (Exceeding 
F  under  HO) 

.27877E-03 


•25564E-03 
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DOLORES  :  MPX  WIND  VS.  TIME 


PREDICTED  VIR  SLTST3  (El  OBS) 

- PREDICTED 

-  OBSERVED 


“•( — - 1 - 1 - 1  t 

24.0  36.0  48.0  60.0  72.0 

T I  ME ( HOURS  ) 


Figure  10.  Predicted  and  observed  maximum  wind  speeds 
for  Dolores  using  three  Independent  variables. 


.0 


12.0 


DOLORES  :  MfiX  WIND  VS.  TIME 


Figure  11.  Predicted  and  observed  maximum  wind  speeds 
for  Dolores  using  four  Independent  variables. 


The  regression  equation  for  Enrique  using  three  independent 
variables  is  of  the  same  form  as  for  Dolores,  except  the  regres¬ 
sion  coefficients  have  the  following  values:  BO  *  968.1, 

B1  ■  -4.46,  B2  ■  0.8934  and  B3  ■  -2.354.  The  four  variable 
equation  has  coefficients  with  the  following  values:  BO  *  277.2, 
B1  -  -1.397,  B2  -  3.583,  B3  »  3.038  and  B4  -  -0.8751.  Plots  of 
the  predicted  vs.  observed  maximum  wind  speeds  via  the  three  and 
four  variable  equations  for  Enrique  are  shown  in  Figures  12  and 
13. 

The  equations  developed  to  predict  rate  of  intensity  changes 
(RIC)  required  the  Inclusion  of  four  variables  to  explain  more 
than  50%  of  the  variance  in  the  case  of  Dolores,  while  there 
wasn't  a  significant  difference  In  the  percent  of  variance 
explained  when  either  three  or  four  variables  were  used  in  the 
case  of  Enrique.  The  analysis  of  variance  tables  for  both  the 
three  and  four  variable  cases  are  shown  in  Table  5  for  Dolores 
and  Table  6  for  Enrique. 

The  regression  equation  for  Dolores  using  three  independent 
variables  is, 

RIC  =  BO  +  B1*IRMAX  +  B2MBMAX  +  B3*ALPHA  (7) 

where  BO  =  -61.29,  B1  *  0.2584,  B2  =  0.0147  and  B3  «  0.029.  The 
four  variable  equation  is, 

RIC  =  BO  +  B1*IRMAX  +  B2*IBMAX  +  B3*IBAV  +  B4*ALPHA  (8) 

where  BO  =  -21.22,  B1  =  0.1694,  B2  =  -0.4194,  B3  =  0.3599  and 
B4  =  0.054. 

The  regression  equations  for  Enrique  are  of  the  same  form  in 
both  cases.  The  coefficients  for  the  three  variable  equation 
are:  BO  =  -90.28,  B1  =  0.34,  B2  =  0.1233  and  B3  =  -0.2275;  while 
for  the  four  variable  equation  they  are  as  follows:  BO  =  -108.6 
B1  =  0.4213,  B2  =  0.1947,  B3  -  -0.0805  and  B4  =  -0.1883.  Plots 
of  the  predicted  vs.  observed  intensity  changes  via  the  three  and 
four  variable  regression  equations  for  Dolores  and  Enrique  are 
hown  in  Figures  14-17. 


20.0  40. 


ENRIQUE  :  MRX  WIND  VS.  TIME 


0.0  12.0  24.0  36.0  48.0  60.0  72.0 

T I  ME ( HOURS ) 


Figure  12.  Predicted  and  observed  maximum  wind  speeds 
for  Enrique  using  three  Independent  variables. 
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ENR I QUE 


:  MRX  WIND  VS.  TIME 

PREDICTED  VIR  SLTST4  (19  OBS) 


24.0 


-I - 1 - 

36.0  48.0 


T I  ME ( HOURS  ) 


60 . 0 
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72.0 


icted  and  observed  maximum  wind  speeds 
sing  four  independent  variables. 
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Table  5.  Analysis  of  variance  table  for  Dolores  (rate  of  intensity  change) 


THREE  INDEPENDENT  VARIABLES 

Source  of 
Variation 

Degrees  of 
Freedom 

Sums  of 
Squares 

Mean 

Squares 

F-Value 

Regression 

. 30000E  01 

.35974E  02 

. 11991E  02 

.51055E  02 

Residual 

.17000E  02 

.39928E  04 

.23487E  01 

Total 

.20000E  02 

. 75901 E  02 

**  GOODNESS  OF  FIT  ESTIMATES  ** 

Percentage  of  Variance  Explained  s  . 
Standard  Deviation  of  Residuals  s  . 
As  a  Percentage  of  Response  Mean  =  . 

47395E  02 
15325E  01 
12513E  04 

FOUR  INDEPENDENT  VARIABLES 

Regression 

. 40000E  01 

.53414E  02 

.13354E  02 

.95012E  01 

Residual 

. 16000E  02 

.22487E  04 

.14055E  01 

Total 

. 20000E  02 

.75901E  02 

**  GOODNESS  OF  FIT  ESTIMATES  ** 


P  (Exceeding 
F  under  HO) 

.10610E-01 


.39351E-03 


Percentage  of  Variance  Explained 
Standard  Deviation  of  Residuals 
As  a  Percentage  of  Response  Mean 


70373E  02 
11855E  01 
96796E  03 


Table  6. 

Source  of 
Variation 

Regression 

Residual 

Total 


Regression 

Residual 

Total 


Analysis  of  variance  table  for  Enrique  (rate  of  intensity  change). 


THREE  INDEPENDENT  VARIABLES 


Degrees  of 
Freedom 

. 30000E  01 

. 15000E  02 

.10000E  02 


Sums  of 
Squares 

.83833E  02 

.29373E  02 

.11321E  03 


Mean 

Squares 

.279744  02 

. 19582E  01 


F-Value 
. 14271E  02 


P  (Exceeding 
F  under  HO) 

.11474E-03 


**  GOODNESS  OF  FIT  ESTIMATES  ** 


Percentage  of  Variance  Explained 
Standard  Deviation  of  Residuals 
As  a  Percentage  of  Response  Mean 


. 74054E  02 
.13993E  01 
. 14394E  03 


. 40000E  01 
.14000E  02 
. 18000E  02 


FOUR  INDEPENDENT  VARIABLES 

.84372E  02  .21093E  02  .10242E  02  .43231E-03 

.28834E  02  .20596E  01 

•11321E  03 


**  GOODNESS  OF  FIT  ESTIMATES  ** 

Percentage  of  Variance  Explained  =  .74530E  02 
Standard  Deviation  of  Residuals  =  . 14351E  01 
As  a  Percentage  of  Response  Mean  =  .14762E  03 
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DOLORES  :  INTENSITY  CHANGE  VS.  TIME 


PREDICTED  VIR  SLTST3  (21  OBS) 
-  PREDICTED 


cs 

•- 

» 


rT- 

0.0 


“I - 

12.0 


24.0 


r 


T 


36.0  4B.0 

TIMECHOURS ) 


60.0 


72.0 


Figure  14.  Predicted  and  observed  rates  of  Intensity 
change  for  Oolores  using  three  independent  variables. 
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DOLORES  :  INTENSITY  CHANGE  VS.  TIME 
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ENRIQUE  :  INTENSITY  CHANGE  VS.  TIME 
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Figure  16.  Predicted 
change  for  Enrique 


and  observed  rates  of  intensity 
using  three  Independent  variables. 
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ENRIQUE  :  INTENSITY  CHANGE  VS.  TIME 
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Figure  17.  Predicted 
change  for  Enrique 


and  observed  rates  of  Intensity 
using  four  independent  variables. 
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4.  BIAS  IN  RESULTS 


The  maximum  wind  speed,  and  the  associated  intensity  changes 

derived  from  them,  used  as  the  dependent  variables  for  the 

regression  equations  are  derived  via  the  Dvorak  technique  and  are 

only  estimates  of  the  maximum  wind  speeds.  The  accuracy  of  the 

Dvorak  technique,  based  upon  independent  evaluations  (Shewchuk 

and  Weir,  1980),  show  that  the  average  error  of  the  satellite 

analysts  was  approximately  0.5  T-number  which  translates  into  an 

estimated  average  maximum  wind  speed  error  of  5-15  kt,  depending 

on  the  position  on  the  Cl  scale.  Results  from  the  Joint  Typhoon 

Warning  Center  (JTWC),  based  upon  post-analysis  intensity  data, 

indicate  a  root-mean-square  (RMS)  error  of  12  kt  as  documented  in 

the  1st  Weather  Wing  Pamphlet  105-10  (1974). 

An  error  analysis  of  the  SLT  predicted  maximum  wind  speeds 

and  intensity  changes  using  the  Dvorak  estimates  as  the  observed 

(ground  truth)  data  are  shown  in  Table  7  for  Dolores  and  Table  8 

for  Enrique. 

I 

I 


4 
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Table  7.  Error  analysis  for  Dolores'  four  regression  cases. 

3  INDEPENDENT  VARIABLES  4  INDEPENDENT  VARIABLES 


0BSVD 

PRED. 

|  error) 

OBSVD 

PRED. 

|  error) 

VMAX 

VMAX 

VMAX 

VMAX 

44.0 

44.7 

.7 

44.0 

44.5 

.5 

48.0 

59.3 

11.3 

48.0 

59.6 

11.6 

66.0 

67.5 

1.5 

66.0 

68.0 

2.0 

66.0 

76.2 

10.2 

66.0 

76.7 

10.7 

66.0 

70.0 

4.0 

66.0 

70.0 

4.0 

68.0 

70.5 

2.5 

68.0 

69.8 

1.8 

70.0 

74.0 

4.0 

70.0 

72.9 

2.9 

80.0 

75.9 

4.1 

80.0 

74.8 

5.2 

94.0 

93.3 

.7 

94.0 

92.0 

2.0 

98.0 

88.5 

9.5 

98.0 

87.7 

10.3 

100.0 

92.6 

7.4 

100.0 

93.0 

7.0 

102.0 

77.6 

24.4 

102.0 

78.7 

23.3 

104.0 

88.9 

15.1 

104.0 

89.7 

14.3 

98.0 

88.6 

9.4 

98.0 

89.0 

9.0 

96.0 

109.1 

13.1 

96.0 

108.9 

12.9 

92.0 

99.9 

7.9 

92.0 

99.9 

7.9 

90.0 

94.8 

4.8 

90.0 

95.0 

5.0 

88.0 

77.7 

10.3 

88.0 

77.5 

10.5 

84.0 

82.5 

1.5 

84.0 

82.1 

1.9 

76.0 

84.3 

6.3 

76.0 

84.4 

8.4 

70.0 

84.3 

14.3 

70.0 

85.6 

15.6 

Sum  of  Squared  Errors 
Mean  Squared  Error 
Root  Mean  Squared  Error 


1985.0290  Sum  of  Squared  Errors  *  1975.2220 

94.5252  Mean  Squared  Error  *  94.0582 

9.7224  Root  Mean  Squared  Error  *  9.6984 


OBSVD 

PRED. 

|  ERROR  | 

RIC 

RIC 

4.0 

2.0 

2.0 

1.7 

2.4 

.7 

1.1 

2.3 

1.2 

.0 

2.2 

2.2 

.0 

1.2 

1.2 

.6 

.9 

.3 

.8 

.3 

.5 

2.1 

.6 

1.5 

2.9 

-1.2 

3.1 

1.7 

.2 

1.5 

.4 

.1 

.3 

.8 

.4 

.5 

.4 

.1 

.3 

-1.3 

-  .3 

1.0 

-  .8 

-  .8 

.0 

-1.7 

-1.1 

.6 

-  .8 

-1.2 

.3 

-  .4 

-1.2 

.8 

-1.7 

-1.6 

.1 

-2.2 

-1.9 

.3 

-5.0 

-1.7 

3.3 

OBSVD 

PRED. 

|  ERROR  | 

RIC 

RIC 

4.0 

2.2 

1.8 

1.7 

2.0 

.3 

1.1 

1.6 

.5 

.0 

1.5 

1.5 

.0 

1.3 

1.3 

.6 

1.7 

1.2 

.8 

1.7 

.9 

2.1 

2.1 

.0 

2.9 

1.6 

1.4 

1.7 

1.3 

.4 

.4 

-  .5 

1.0 

.8 

-1.1 

1.9 

.4 

-1.0 

1.4 

-1.3 

-  .8 

.4 

-  .8 

-  .5 

.3 

-1.7 

-1.1 

.6 

-  .8 

-1.4 

.6 

-  .4 

-1.0 

.5 

-1.7 

-1.1 

.5 

-2.2 

-2.0 

.2 

-5.0 

-3.4 

1.6 

Sum  of  Squared  Errors  s  39.9289 
Mean  Squared  Error  =  1.9014 
Root  Mean  Squared  Error  =  1.3789 


Sum  of  Squared  Errors  s  22.4912 
Mean  Squared  Error  =  1.0710 
Root  Mean  Squared  Error  =  1.0349 
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Table  8.  Error  analysis  for  Enrique's  four  regression  cases. 


3 

INDEPENDENT  VARIABLES 

4  INDEPENDENT  VARIABLES 

OBSVO 

PRED. 

OBSVD 

PRED. 

VMAX 

VMAX 

1  ERROR  | 

VMAX 

VMAX 

|  ERROR  | 

56.0 

72.4 

16.4 

56.0 

72.5 

16.5 

58.0 

73.7 

15.7 

50.0 

67.7 

9.7 

60.0 

73.0 

13.0 

60.0 

66.4 

6.4 

62.0 

67.2 

5.2 

62.0 

67.4 

5.4 

66.0 

82.5 

16.5 

66.0 

70.0 

4.0 

78.0 

87.6 

9.6 

78.0 

73.7 

4.3 

83.0 

83.9 

.9 

83.0 

76.0 

7.0 

96.0 

90.9 

5.1 

96.0 

95.0 

1.0 

105.0 

99.3 

5.7 

105.0 

112.7 

7.7 

112.0 

123.3 

11.3 

112.0 

122.4 

10.4 

124.0 

119.3 

4.7 

124.0 

112.9 

11.1 

125.0 

117.6 

7.4 

125.0 

109.6 

15.4 

125.0 

99.6 

25.4 

125.0 

103.3 

21.7 

125.0 

105.5 

19.5 

125.0 

110.2 

14.0 

125.0 

103.2 

21.8 

125.0 

100.6 

24.4 

125.0 

123.7 

1.3 

125.0 

121.7 

3.3 

124.0 

136.2 

12.2 

124.0 

128.4 

4.4 

117.0 

136.1 

19.1 

117.0 

131.5 

14.5 

104.0 

131.0 

27.0 

104.0 

128.1 

24.1 

Sum  of  Squared  Errors 

s 

4087.2030 

Sum  of  Squared  Errors 

= 

3154.8870 

Mean  Squared  Error 

= 

215.1160 

Mean  Squared  Error 

= 

166.0467 

Root  Mean 

Squared  Error 

S 

14.6668 

Root  Mean  Squared  Error 

= 

12.8859 

OBSVD 

PRED. 

OBSVD 

PRED. 

1  1 

RIC 

RIC 

|  ERROR  | 

RIC 

RIC 

|  ERROR  | 

.8 

1.9 

1.1 

.8 

2.0 

1.1 

.8 

1.2 

.3 

.8 

1.1 

.3 

.8 

1.3 

.4 

.8 

1.2 

.4 

.8 

1.0 

1.0 

.8 

1.9 

1.0 

1.7 

2.6 

1.0 

1.7 

2.4 

.7 

5.0 

3.2 

1.8 

5.0 

2.9 

2.1 

2.1 

2.8 

.7 

2.1 

2.7 

.6 

5.4 

2.7 

2.7 

5.4 

2.9 

2.5 

3.8 

3.4 

.3 

3.8 

3.9 

.1 

2.9 

4.2 

1.3 

2.9 

4.2 

1.3 

2.5 

2.8 

.3 

2.5 

2.7 

.2 

.4 

1.3 

.9 

.4 

1.1 

.7 

.0 

-  .5 

.5 

.0 

-  .4 

.4 

.0 

-  .8 

.8 

.0 

-  .6 

.6 

.0 

-1.4 

1.4 

.0 

-1.4 

1.4 

.0 

-  .8 

.8 

.0 

-  .8 

.8 

-  .3 

-1.8 

1.6 

-  .3 

-2.0 

1.7 

-2.9 

-2.4 

.5 

-2.9 

-2.4 

.5 

-5.4 

-2.9 

2.5 

-5.4 

-2.9 

2.5 

Sum  of  Squared  Errors 

s 

29.3720 

Sum  of  Squared 

Errors 

= 

28.8347 

Mean  Squared  Error 

= 

1.5459 

Mean  Squared  Error 

= 

1.5176 

Root  Mean  Squared  Error 

= 

1.2434 

Root  Mean  Squared  Error 

= 

1.2319 
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5.  CONCLUSIONS  AND  RECOMMENDATIONS 


The  foregoing  predictive  equations  are  not  meant  to  be  used 
as  prognostic  equations,  rather  they  are  a  measure  of  the 
feasibility  of  the  SLT's  usefulness  as  a  forecast  tool.  The 
significant  differences  for  the  two  storm's  set  of  regression 
coefficients  leads  to  the  question  as  to  which  set  of  coeffi¬ 
cients  to  use  for  a  given  storm.  Gentry  et  al.  (1978)  found  it 
necessary  to  stratify  the  storms  based  upon  their  current  inten¬ 
sity  estimate,  with  an  equation  for  weak  storms  (Vmax<.65  kt) 
and  one  for  more  intense  storms.  The  stratification  of  the 
storms  is  useful  only  for  intensity  forecasts.  More  work  needs 
to  be  done  with  independent  data  sets  before  any  statistically 
significant  results  may  be  ascertained. 

In  the  case  of  maximum  wind  speed  estimation,  the  storm  to 
storm  variability  in  the  amount  and  structure  of  the  convective 
activity  appears  to  be  too  large  for  the  development  of  an 
objective  technique  utilizing  SLT-derived  parameters  to 
accomplish  such  a  task.  This  result  agrees  with  that  of  Brody 
and  Tsui  (1982)  in  their  attempt  to  predict  maximum  wind  speeds 
using  the  technique  described  by  Gentry  et  al.  (1978). 

The  outlook  for  the  development  of  an  objective  technique 
for  predicting  intensity  changes,  however,  appears  promising. 
SLT-derived  parameters  used  in  conjunction  with  parameters 
derived  from  the  tropical  cyclone  intensity  forecast  program 
(Brody  and  Tsui)  using  the  algorithm  developed  by  Gentry  et  al. 
may  provide  the  answer.  Other  parameters  such  as  latitude  of  the 
storm,  time  of  year,  and  oceanic  basin  ^ay  help  in  explaining 
more  of  the  observed  variance.  Inclusion  of  the  current 
intensity  estimate  into  the  forecast  scheme  may  result  in  a  more 
accurate  intensity  change  estimate,  although  with  somewhat  less 
objectivity. 
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Smith  (1972)  used  low  and  high  level  cloud  motions  derived 
from  ATS-11  satellite  images  to  construct  fields  of  radial  and 
tangential  wind  components,  areal  relative  vorticity,  and  areal 
divergence  for  various  times  during  the  life  cycle  of  Hurricane 
Celia  (1970).  Using  his  values  for  the  low  cloud  level  and  the 

relations  described  by  Lahlrl,  (l.e.,  v*V/  *  -a  and  Vr/Vo  “  _a) 
the  computed  ranged  for  the  spiral  constants  for  Hurricane  Celia 
were  found  to  be  0.2  <  a  <  0.5  for  large  radii.  This  translates 
to  Inflow  angles,  a,  ranging  from  11.30  <  a  <  26.60  for  the  wind 
field  spirals.  Comparing  this  with  the  range  of  SLT-derived 
cloud  band  spiral  constants  and  inflow  angles  of  0.3  <  a  <  0.8 
and  16.7o  <  a  <  38.70  respectively.  It  appears  that  the  spirals 
traced  out  by  the  cloud  bands  have  a  much  larger  inflow  angle 
than  that  of  the  wind  field.  The  assumption  made  by  Lahlrl,  that 
the  cloud  band  spirals  can  be  used  as  a  tracer  of  the  wind  field 
spirals  appears  to  be  invalid. 

The  current  SLT  has  only  limited  effectiveness  in  isolating 
spiral  banding  features  from  the  large  and  sometimes  amorphous 
cloud  masses  associated  with  tropical  cyclones.  The  current 
method  isolates  a  spiral  band  by  enhancing  the  satellite  image. 
The  spiral  band  must  be  clearly  distinguishable  from  cloud 
features  not  associated  with  the  spiral  band.  The  resulting 
transformed  spiral  band  must  then  be  readily  subsetted  (excluded) 
from  the  cloud  masses  not  associated  with  the  spiral  band  so  that 
the  least  squares  line  will  be  applied  to  the  band  only.  This 
process,  at  times,  cannot  isolate  obvious  spiral  features  neces¬ 
sary  to  derive  the  SLT  parameters. 

A  new  SLT  is  being  developed  which  will  allow  the  user  to 
define  a  spiral  feature  using  the  graph  pen  and  pad.  This 
eliminates  the  necessity  to  clearly  isolate  the  spiral  band,  as 
with  the  current  SLT,  and  does  not  require  the  image  to  be 
transformed,  pixel  by  pixel,  into  spiral  coordinates,  but  rather 
a  spiral  equation  can  be  fitted  to  the  curve  traced  out  with  the 
graph  pen.  This  allows  for  the  same  statistical  parameters  to  be 
calculated  as  well  as  the  objective  determination  of  the  spiral 
constant,  center,  and  offset  radn.  More  important,  however,  is 
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that  the  new  method  allows  for  the  use  of  the  FFT's,  at  specified 
points  both  along  and  normal  to  the  spiral  band,  to  determine 
information  concerning  cloud  size  and  spacing  previously  unavail¬ 
able  due  to  distortions  caused  by  the  current  SLT  transformation 
algorithm . 

The  results  of  this  project  are  sufficiently  promising  to 
warrant  further  investigation  into  the  development  of  an 
objective  technique  to  determine  intensity  changes  employing 
parameters  derived  from  the  new  SLT  and  those  derived  from  the 
tropical  cyclone  intensity  forecast  algorithm. 
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